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During the last years, the ALPHA collaboration has been developing and implementing a method 
based on Heavy Quark Effective Theory (HQET) to compute B-mesons observables through lat- 
tice simulations. Thanks to a non-perturbative matching to QCD, the theory is renormalizable at 
any order of the heavy quark mass expansion. In order to extract precisely the relevant matrix 
elements and masses, we use all-to-all propagators and solve an generalized eigenvalue problem 
(GEVP). We have shown in the quenched approximation that quantities like the i>-quark mass mb, 
the heavy-light decay constant(s) or the Z?-meson spectrum can be computed precisely beyond the 
static approximation (including the first corrections in 1 /m\j). More recently, we have started to 
include the sea quark effects, by working with Nf = 2 light flavors of dynamical fermions [1] . 
The computation of the matching parameters is almost finished, but concerning the extraction of 
the hadronic quantities for which we use some CLS ensembles [2], only one lattice spacing has 
been analyzed so far. In this proceeding we report on the status of this project and present some 
preliminary results. The strategy is sketched in the following figure; 
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1. Introduction 

The search for new physics is currently limited by the size of the theoretical uncertainties, 
mainly due to the errors affecting the hadronic quantities. This is specially true in the B physics 
area, where in many cases the errors are largely dominated by quantities like heavy-light decay 
constants or bag parameters, which are obtained by lattice simulation. If the recent progress of 
the lattice community in the light quark sector are very impressive, lattice simulations around the 
Zj-quark mass are still difficult. Various strategies have been developed (see e.g. [3-5]) based on 
effective theories (NRQCD, HQET), on relativistic formulation (Fermilab, RHQ) or a combination 
of both. As already mentioned in the abstract, the strategy followed by the ALPHA collaboration 
is based on non-perturbative HQET in the static approximation and at the 1/mb order (see [6] for a 
pedagogical introduction). This is a theoretically sound framework, in which the systematic errors 
are well under control. Based on our experience in the quenched approximation, we believe that 
the total uncertainty affecting the results can be made small enough to have an impact in the search 
for new physics, and maybe unveil BSM effects. This method has been tested in the quenched 
approximation to study the meson spectrum [7], to determine the ^7-quark mass [8] and the 
decay constant /b^ [9]. The computation is essentially done in two steps: as explained in section 2, 
the HQET parameters are obtained from a non-perturbative matching to QCD in a small volume, 
while the hadronic matrix element and masses are obtained from a subset of ensembles generated 
within the CLS effort (see section 3). We present our conclusions together with our prehminary 
results in section 4. 

2. Computation of the relevant HQET parameters 

At the 1/mb order, we write the HQET Lagrangian in the following way: 

^HQEt(jc) = .^%xai{x) - G)idn^kin(jc) - ©spin ^spin (jc) , (2. 1) 

^statW = FhW^O V^hW , ^kinW = VAh(x)DVhW , ^spinW = Wh{x)<^ '^Whix) ■ (2.2) 

We are also interested in the time component of the heavy-light axial current Aq. Considering only 
the terms which contribute to zero spatial-momentum correlation functions, we write it as 

^^(x) = Z«Q^^ [Ar (x) + c2Ui^)(x)] , (2.3) 
W = FiW ^75T^-(Vf- Vf ) VAh(x) , Af\x) = ?i W7075 V^h (x) , (2.4) 

where Vf denotes the symmetric derivative. The computation of (Okinj^spim^A^^^jCyiJ^ and mbare 
is done following the strategy presented in [1, 10]. In a volume of linear space extent L\ ~ 0.5 
fm, where the Z7-quark can be simulated with discretization effects under control, we compute a 
set of observables <I>;=i 5 at four different lattice spacings and extrapolate them to the continuum 
^ In this set of simulations called 5i in the figure, the light quark masses are set to 0, while for 

'in the continuum and large volume limits, $1 is proportional to the meson mass and $2 to the logarithm of the 
decay constant, respectively. 4>3 is used to determine the counter-term of the axial current, and 4>4 5 for the determination 
of the kinetic and magnetic term, respectively. 
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the (RGI) heavy quark mass M we have chosen nine different values covering a range between the 
charm to the bottom mass. In another set of simulations called S2, where we use the same value of 
the physical volume, we compute the corresponding quantities in the effective theory (at the 1 /mb 
order) and match them to their QCD counterpart. This matching can be written in the following 
way: 

^f^^{Li,M,0) = r]i{Li,a)+'£(pij{Li,a)o)j{M,a) , (2.5) 

where T] and (p are computed by lattice simulations for different values of the lattice spacing a. In 
other words the matching equations determine the set of parameters a) = (p^^ [<J>QCD — rj]. We then 
compute rj and (p in a larger volume of space extent L2 = 2Li, and using the parameters a){M,a) 
determined in the previous step we compute the observables <I>(L2,M,0) according to a formula 
similar to eq. (2.5) but where Li is replaced by L2 and w by O) (note that in this step called S3 
the continuum limit can be taken because the divergences cancel out exactly). This procedure can 
then be re-iterated until the volume reached is large enough for finite size effects to be negligible, 
typically around (2fm)^. In practice, it turns out that three different volumes are enough (Li ,L2 and 
the large volume one). Thus, the HQET parameters that can be used in large volume simulations 
(denoted as 54) are given by (o{M,a) = (p^\L2,a)[(i>{L2,M,0) - ri{L2,a)] . 

3. Extraction of HQET hadronic matrix elements and preliminary results 

In the large volume limit, at the first order of the 1 / nib expansion our main observables are 
given by 

'MB = mbare + E''"' + (Okin^"^" + £Ospi„£^P'" , (3. 1) 

mB-mB* = ^Wspin^'P'", (3.2) 
log(a3/2/BvW2) = log(Z^Q^T)+log(a3/2p^'-') + Z7f' amq 

+(Ok,n;^'^" + Wspin/^^P'" + ci^'/^^" , (3.3) 

where b^^^ is an improvement coefficient. The HQET energies £'stat ^-km ^■spin matrix el- 

ements p^'^^jP^''* ,p'''",p*'P'" have been measured on a subset of configuration ensembles produced 
within the CLS effort [2] with Nf = 2 flavors of (a) -improved Wilson-Clover fermions. We have 
employed a GEVP analysis in order to control the excited states contamination. The chiral extrap- 
olations are performed using pion masses down to 250 MeV (and the coupling gB*B;t computed 
in [11]). In principle these quantities should be calculated at various values of the lattice spacing in 
order to take the continuum limit. This continuum extrapolation will be left for future work, instead 
here we just compute the observables at one (quite small) value of the lattice spacing a = 0.07 fm. 
From our experience in the quenched approximation, we do not expect the discretization effects to 
be visible within our present accuracy. More details of this computation can be found in [1]. 

4. Preliminary results and conclusion 

For the b-quark mass, including the l/wib terms, we find : 

. -f. ^ -„ HQET 

m^^im^^) =4.276(25),„(50),tat+renorm(?)« GeV. (4.1) 

Nf— 2 
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where the first error comes from the uncertainty on ro, while the second error includes the statis- 
tical error on aE^^^, the uncertainty on the chiral extrapolation and the error on the quark mass 
renormalization constant. At the same order our result for the heavy-light decay constant /b reads 

/bn?=2 = 178(16)(?)aMeV, (4.2) 

where the first error includes the statistical uncertainty on matrix elements, the systematics coming 
from chiral extrapolation and the uncertainty on the physical scale of tq. In both cases the "(?)a" 
indicates that a continuum limit is not yet performed, but we expect the discretization effects to be 
small. These first results are promising and once we have controlled cut-off effects by simulation 
at several lattice spacings, we plan to extent our project to the computation of other hadronic quan- 
tities like the B — B mixing or B TV semileptonic form factors as well as the spectrum of hadrons 
with a b-flavor. 

I am indebted to my colleagues of the ALPHA collaboration, in particular I would like to 
thank F. Bemardoni, B. Blossier, J. Bulava, M. Delia Morte, M. Donnellan, P. Fritzsch, J. Heitger, 
G. von Hippel, T. Mendes, H. Simma and R. Sommer for many stimulating discussions and contri- 
butions to this work. 
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